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Abstract. Given a non-uniform criss-cross partition of a rectangular domain O, we analyse 
the error between a function f defined on f2 and two types of C 1 -quadratic spline quasi- 
interpolants (QIs) obtained as linear combinations of B- splines with discrete functionals as 
coefficients. The main novelties are the facts that supports of B-splines are contained in £1 
and that data sites also lie inside or on the boundary of £1. Moreover, the infinity norms 
of these QIs are small and do not depend on the triangulation: as the two QIs are exact 
on quadratic polynomials, they give the optimal approximation order for smooth functions. 
Our analysis is done for f and its partial derivatives of the first and second orders and a 
particular effort has been made in order to give the best possible error bounds in terms of 
the smoothness of f and of the mesh ratios of the triangulation. 
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1 Introduction 

Given a non-uniform criss-cross partition of a rectangular domain f2, we analyse the error 
between a function / denned on Q and two C 1 quadratic spline quasi-interpolants (abbr. 
QIs), denoted 52 and W|, obtained as linear combinations of B-splines with discrete co- 
efficient functionals. The first operator S2 was described by the second author in |llj|14j 
and the second one is a slight modification of the operator W2 introduced by Chui and 
Wang in 4_, and also studied by Chui and He in j2. , Wang and Lu in JH] and by the first 
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author in jS], [7]. 

With respect to previous papers, we note the following facts : we introduce B-splines with 
supports contained in Q and data sites lying inside or on the boundary of ft, so we do 
not need extra values outside the domain. This can be useful in certain practical problems 
where these data are not available. Moreover, we show that the infinity norms of these QIs 
are small and do not depend on the triangulation. As they are exact on the space P2 of 
quadratic polynomials, it is well known that they give the optimal approximation order for 
smooth functions. 

Another important and very useful property of QIs is that the construction of these opera- 
tors do not need the solution of any system of equations. It is particularly attractive in the 
bivariate case where the number of data sites can be huge in practice. 
Though the QIs do not interpolate / at data sites, it can be observed that errors are quite 
small at that points. Actually, a superconvergence phenomenon can often be observed at 
some specific points. Moreover, the global behaviours of QIs and of their derivatives is quite 
close to those of the function / (see e.g.jH]). 

Our error analysis is done for / and its partial derivatives of the first and second orders 
and a particular effort has been made in order to obtain sharp error bounds in terms of the 
smoothness of / and of the characteristics of the triangulation, in particular local mesh ra- 
tios. Such a program can be developed thanks to the good properties of quadratic B-splines 
described in 13 j. It is true that we do not get the best error constants, which is a rather 
technical task, however, we obtain a reasonable order of magnitude of these constants. This 
can be useful in the practical applications that we want to develop elsewhere. 
Here is an outline of the paper: in Section 2, we recall the main definitions on the B-splines 
on criss-cross triangulations that we use in the definition of quasi-interpolants . In Section 
3, we describe the two quadratic spline QIs. In Section 4, we give error estimates of the 
infinity norms of / — Q, where Q = S2 or W£, when / G C s (f2), with < s < 3. In Section 
5, we give error estimates on first derivatives \\D r,s (f — Q)\\ OQ , r + s = 1, in 0, and on second 
derivatives \\D r ' s (f — <3)l|oo> r + s = 2, inside triangular cells of the triangulation, since Q is 
only C . They are expressed in terms of moduli of smoothness with respect to the length 
h/2, where h is the maximal steplength of the given partition of the domain. 

2 Quadratic B-splines on a bounded rectangle 

In this Section, we first introduce C 1 quadratic B-splines generating the spline space in 
which we approximate functions. Then, in the following section, we will define the two 
quasi-interpolants S2 and W%. 

Let f2 = [a, b] x [c, d] be a rectangle decomposed into mn subrectangles by the two partitions 

X m = < i < m}, Y n = {yj,0 <j< n}, 

respectively of the segments / = [a, b] = [xq, x m ] and J = [c, d] = [yo, y n \. We also introduce 
the double knots x_i = x , y_i = y , x m+1 = x m , y n+ i = y n . 
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The so-called criss-cross triangulation T mn of is defined by drawing the two diagonals in 
each subrectangle Qij = [xi-i,xi\ X \yj-i,Vj]. We need the two following sets of indices: 

JC m n = {(i,j) :0<i<m + l, 0<j<n + l}, 

K-mn = :l<i<m, l<j<n}. 

We set hi = Xi-Xi-i, kj =yj-yj-i, S{ = \{xi-\ + Xi), tj = \{Vi-i + Vj) for (t, j) G /C mn . 
We denote by = (xi,yj), 0<i<m, 0<j<n} the set of vertices of subrectangles 

and by {Mjj = (si,tj), G K, mn } the set of their centers, of midpoints of boundary 

subintervals and of vertices of 0,. 

Let B mn := {-By, (i, j) £ K, mn } be the collection of (m + 2)(n + 2) B-splines generating the 
space S2{T mn ) of all C 1 piecewise quadratic functions on the criss-cross triangulation T mn , 
associated with the partition X m x Y n of the domain O. There are mn B-splines associated 
with the set of indices )C mn , whose restrictions to the boundary T of are equal to zero. 
They were also introduced in jSJj^lSl- To the latter, we add 2m + 2n + 4 boundary B-splines 
whose restrictions to V are univariate quadratic B-splines. Their set of indices is 

ICmn ■= {(*, 0), (t, n + 1), < i < m + 1; (0, j), (m + 1, j), < j < n + 1}. 

The BB (=Bernstein-Bezier)-coefficients of inner B-splines {B^, 2 < i < m — 1, 2 < j < 
n — 1} are given in The other ones can be found in the technical report ^3] an d 

in ^21- The B-splines are positive and form a partition of unity (blending system). The 
boundary B-splines are linearly independent as the univariate ones. But the inner B-splines 
are linearly dependent, the dependence relationship being: 

^ i-iy^hikjBij = o. 

{i,j)EKmn, 

Although S mn is not a basis of 52(^mn)) this fact has no influence on the definition and 
properties of QIs. The support of Bij is denoted by Sy: for inner B-splines, it is a non- 
uniform octagon. The set B mn can also be defined in the following way. Define the extended 
partitions 

and 

Y n = Y n U {y_2> Vn+li Vn+2ii 

where x_ 2 < X-\ < x , x m < x m +i < x m +2,_ y~- 2 < V-i < Vo, Vn < y n +i < Vn+2, 
and the corresponding criss-cross triangulation T mn . We also put ho = xo — x-\, h m+ \ = 
x m+1 - x m , k =y - y_ 1} kn +1 = y n± i - y n . 

We consider the collection B mn := {By,(i,j) £ K, mn } of the (m + 2)(n + 2) "classical" 
B-splines with octagonal supports such that n int(Q) ^ |3]. 

We note that Bij = B^ for inner B-splines. Using the BB-coefficients of both families 
B mn and B mn , one can derive the expressions of the new boundary B-splines in function of 
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"classical" B-splines. For this purpose, we need the following notations, for 2 < i < m and 
2 < j < n : 



hi / hi— i kj / kj—\ 

fj . = = — q- ■ f — - T • = - 

hi-i + hi' 1 hi_i + hi % kj^i +kj' 3 kj-i + kj 



In addition, we need the particular values : 



hi / h m _ k\ , k n 



° X ~ T 717 ' a m+l — 7 ~T ' T l — f 7T~ ' T n+l 



ho + h\ h m + h m +i ko + k\ k n + k n+ \ 

and (Ji = a' m+1 = n = r' n+1 = 1, whence cr^ = a m+ i = t[ = T n+1 = 0, since h = h m+1 = 
ko = K + i = 0. 



The first boundary layer of B-splines along the horizontal edge Aoo^mO is defined by 
-Boo = Bqo, Bio = w^{Bio - =^oo)> B i0 = =^B i0 , 2 < i < m — 1, 



B m o — =7(B m fi - ^V >+1 B m+ i o), B m+ \fi — =i 1 =- B m+ ifl- 

1 m + l m + 1 1 

In the same way we obtain, along the vertical edge AooAo n , 

Boi = ±(Boi ~ J#o,o), Boj = ±-B 0j , 2 < j < n - 1, 

Bon = =-(Bon — Sr^i^Crt+l ), Bo^n+l = ="=? #0,n+l- 

Similar formulas hold for boundary B-splines along the edges A m0 A mn and Ao n A mn : 

{B ijn+1 ,0<i<m+l} and {B m+1J , < j < n + 1}. 

The restrictions of all these B-splines to the boundary of Q are classical univariate quadratic 
B-splines. The second boundary layer of B-splines along the horizontal edge Aoo^mO is 
defined by 

Bn = B n - ^-Bi,o - ^#oi + flfT^oo, 
and similar formulas for B m \ , B\ n and B mn . 

Finally we define the second layer of B-splines along the vertical edge AqoAou, 

— T ' — — a' — 

Bn = B a - ^-B ifl , 2<i<m-l, B ±j = B ±j - ^-B ,j, 2 < j < n - 1, 

Tl CTl 

and similar formulas for the collections: 

{B in , 2<i<m-l} and {B mj , 2<j<n-\}, 
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with the ratios and instead of ^ and ^ respectively, in the formula defining 

T n+l a m+l T l ff l 

B n . 

Remark : note that many coefficients can be simplified, for example 



'i 



i T 'i _ ko T n+l _ k n+ i 



cri hi p' m +i h m ti ki r' n+1 k n 

Error analyses given in sections 3 and 4 below are based on the Bernstein Bezier represen- 
tation of B-splines on the triangulation. The associated techniques are described e.g. in [Q, 

El, cm. 

3 Quasi-Interpolants exact on P2 

We now define the two quadratic spline quasi-interpolants S2 and that we want to 
study. Moreover, we give uniform bounds on their infinity norms. 

3.1 The quasi-interpolant S 2 

For the definition of S2, we need the notations Oi and r,- given above in Section 2. Then we 
define: 

° 2 i<+l °i{< + l? _ _ T J T J+1 , _ T 3( T 'j+l? 



bij = 1 — (flij + Cj + CLj + Cj), 

with a = c = a m+ i = c m+ i = a = c = a n+1 = c n+ i = and b = b = b m+ i = b n+ i = 1. 

The data sites for S 2 are the (m + 2)(n + 2) points of the set 

V mn := {Mij = (si,tj), £ K m n} ■ 
The quadratic spline quasi-interpolants £2 jl2j |14j is defined as follows: 

m+l n+ 1 
i=0 j=0 

with coefficient functionals given by 

(1) mj{f) = hjfiMij) + ai f(M^i tj ) + Cif(M i+1J ) + S j /(M ij -_i) + c j f(M i>j+1 ). 

It is exact on P2 and its infinity norm is uniformly bounded independently of the triangu- 
lation T mn of the domain. Indeed, since 

(2) |aj|, \ci\, \a,j\, \cj\ < 1/2 and \bij\ < 3, 
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then it is clear that 

Halloo < 5. 

We notice that the number of data sites requested by S2 is equal to 

(3) N s = ran + 2m + 2n + 4. 

3.2 The quasi-interpolant W 2 * 

The second quasi-interpolant W£ here analysed is a modification of the QI W2 derived by 
Chui-Wang |I] , which is also exact for P2 . The latter is defined in terms of classical B-splines 
{Bij} on the triangulation T mn . Given the values of a function / at the (m + 3)(n + 3) 
points Aij = (xi,yj) , — 1 < i < m + l, —1 < j < n + 1 (among which those having one 
extra abscissa or ordinate are outside Q) and the (m + 2)(n + 2) points Mij, intersections of 
the diagonals in the subrectangles with vertices Ai—ij—i,Aij,Ai—ij,Ai j—i (among which 
a number also lay outside the domain), the Chui-Wang QI is defined by: 

m+l n+1 

i=o i=o 

with coefficient functionals defined by 

^■(/) = 2/(M M ) - ^{/(Ai-u-!) + /(^-ij) + /(^ij-i) + /(Ai)]- 

In that case, the number of data sites is equal to 

N\y = 2mn + 3m + 3n + 9. 

Now if we set x_ 2 = sc_i = x , x m+2 = x m+1 = x m , y_ 2 = y~i = yo, y n +2 = y n +i = y n 
and if we use the B-splines defined in Section 2, with supports Ejj included in the 
domain f2, we can define the modified Chui-Wang QI as follows: 

■m+l n+1 
j=0 j=0 

where the coefficient functionals are: 

(4) M * (/) = 2/(Af&) - ^[/(4*-i,,-i) + f(AU tj ) + f(Al^) + /(A*j)] , 
with the new data points : 





A* 


- A- 


for < i < m, < j < n, 


(5) 


A* 


= ^,0) 


A i,n+1 = A i,n, -1 < i < m + 1, 




A* 


= A j , 


A m+ i j = A m j, —1 < i < m + 1, 




M* ■ 


= M i]:i 


for < i < m + 1, < j < n + 1 
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The number of data sites requested by is equal to 

(6) Nw = 2mn + m + n + 1, 

and they all lie inside the domain Q or on its boundary. 
From \n*j(f)\ < 3| |/| |oo, we can immediately deduce: 

1 1^2 I loo < 3 

for all non-uniform triangulations T mn of the domain 0. 

We remark that both S2 and W£ are local schemes, because for (x, y) G f2, the values 
S2f(x,y) and W%{x,y) only depend on those of / in a neighbourhood of (x,y). 
If < r < m — 1 and < s < n — 1 are integers such that x G [x r , x r +i], y G [y s , then 
(x, y) will belong to one of the four triangular cells Tg of T mn , £ = 1,2,3,4, labelled as in 
Fig. □ 

Each triangle Tg, i = 1,2,3,4, is covered by exactly seven supports of B-splines Sy. In 
Table 1 below, we report the set K(Tg) of indices of such B-splines, as functions of r and s, 
i.e. K(T e ) = E i:j n int(T e ) ± 0}. 

Therefore, if (x, y) G Tg, then: 

(i,j)£K(T e ) 

Wif(x,y)= ^U)Bij{x,y). 



y s +i 

Vs 



Figure 1: Four different kinds of cells in T n 
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Ti 


T 2 


T 3 


T 4 




r, s — 1 


r — 1, s — 1 


r — 1, s — 1 


r, s — 1 




r — 1, s 


r, s — 1 


r, s — 1 


r + 1, s — 1 




r, s 


r — 1, s 


r + 1, s — 1 


r — 1, s 




r + 1, s 


r, s 


r — 1, s 


r, s 




r-1, s + 1 


r + 1, s 


r, s 


r + 1, s 




r, s + 1 


r-1, s + 1 


r + 1, s 


r, s + 1 




r + 1, s + 1 


r, s + 1 


r, s+1 


r + 1, s + 1 



Table 1. 



4 Error analysis for functions 

In this section we analyse the errors / — £2/ and / — W£f for / G C S (Q), < s < 3. We 
need to introduce the following notations: 

h = max{/ij, kj}, 5 = min{/ij, kj}; \ \ • ||oo,fi = || • ||o, = supremum norm over f2; 

= ^(aLaa) = ^ with | a | = ^ + ^ w(£)8/j t) = max{uj (D a f, ^ | a | = s } ; 

||(x,y)|| = (x 2 + 1/ 2 ) 1 / 2 ; e a (x,y) = x ai y a2 = monomial of total degree \a\, 
where the modulus of continuity of ip £ C(Q) is given by: 

u{ip,t) =max{|^(M) - i/j(P)\; M,Pe Sl,\\MP\\ < t}. 
We denote by Q the generic quasi-interpolant defined by: 

m+l n+1 

(7) q/=EEm/)% 

1=0 jr'=0 

where the coefficient functionals are defined by Xij = Hij (1) when Q = S2, and \j = //*■ 
(4) when Q = W 2 *. 

Theorem 1. (Error bounds for continuous functions). There exists a constant Cq > 0, 
with C < 20.5 for Q = S 2 and C < 12 for Q = W% , such that, for f e (7(0) 

||/-Q/||n<C w(/,^). 
Proof. We consider some closed triangular cell T of T mn , for which 

II/-Q/I|q = II/-Q/I|t. 
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T is one of the four triangles depicted in Fig. ^ For the sake of simplicity we assume that 
T = T3. Since Q reproduces P2, for any P £ T, we can write f(P) = Yl(ij)eK(T) f(P)Pij- 
For Q = S 2 , we can write the following inequality: 

\(S 2 f-f)(P)\ < ]T B ij {\b ij \\f(M ij )-f(P)\ + \a i \\f(M i ^ j )-f(P)\ + 

(i,j)eK(T 3 ) 

IdWfiMi+u) - f(P)\ + faWfiMij-t) - f(P)\ + \c,\\f(M l , J+1 ) - f(P)\} 

Assuming that the origin lies at the midpoint of the lower edge of T3, then this triangle 
can be decomposed into two equal subtriangles by the y-axis. By the symmetry of the 
problem, it is sufficient to consider the case when the point P = (x, y) lies in the right 
triangle. Therefore the coordinates satisfy < x + y < ^. We shall now use the following 
simplified notations: there are seven B-splines whose supports intersect int(T) and we 
denote their centres by {M k , 1 < k < 7}, with Mi = M r>s+ x, M 2 = M r _i )S , M 3 = M r , s , 
M 4 = M r+ i iS , M 5 = M r _i iS _i, M 6 = M rjS _i, M 7 = M r+ i iS _i. Each central point M k has 
four neighbours N k , S k , E k , W k (for North, South, East and West positions) involved in the 
coefficient functional [i k . The biggest constants being obtained for k = 1,2, 5, we only detail 
one of these cases, for example k = 1. Then, we obtain the following upper bounds for the 
various distances involved in the majoration : 

||PM 1 ||<VIo|<4^ \\P Nl \\<V26^<6^ \\P Sl \\<V2^<2^ 

11 hi - 2 _ 2> 11 111 _ 2 ~ 2 

Using inequalities ©, we see that the coefficient of the B-spline B>\ = B TiS+ i, whose support 
is centered at M\ is first bounded above by 

3u(f, llPMxH) + 1 (u,(/, \\PNiW) +u(f, \\PS 1 \\)+u(f, \\PE4) +u(f, ||PWi||)) , 
then, using the above upper bounds on distances, we see that it is bounded above by : 

[12 + 1(6 + 2 + 4 + 5)] cj(f, ±) = 20.5o)(/, \). 
Finally, since Ylu j)eK(T 3 ) Bij = 1> we obtain, for all P E T : 

I (S 2 f-f)(P)\<20.5u(f,~). 
which proves that ||/ - S 2 f\\ < 20.5 §). 
Similarly, for Q = WJ, we can write the following inequality: 



(i,i)6if(T 3 ) 

- /(P)| + |/(^+i )i+1 ) - /(P)|)}. 
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We now compute upper bounds for the distances involved in the case when = (r, s+1). 
The central point M x = M rjS+ i has four neighbours N W 1 , NE X , SW 1 , SWi (for North- West, 
North-East, South- West and South-East positions) involved in the coefficient functional /x*, 
thus we obtain : 

||PMi||<VlO^<4^ ||PWi||<2V5^<5^, \\PNE4 < Wf^ < 5^, 

||PSWi||<2v^^<3^, \\PSW 1 \\<V5 f ^<3^. 

We see that the absolute value of the coefficient of the B-spline B\ = B r)S+ \ whose support 
is centered at Mi is first bounded above by 

2oj(f, \\PM4) + (u(f, IIPWill) + u(f, \\PNE 1 \\)+u(f, \\PSW4) +u{f, \\PSE 1 \\)) , 
then, using the above upper bounds on distances, we obtain as upper bound : 

[8 + ±(5 + 5 + 3 + 3)M/,£) = 12 W (/,£), 
and, finally we obtain, for all P G T : 

\Wf-f){P)\<\2 U <J±) 
which proves that ||/ - W£f\\ < 12 §). ■ 

Theorem 2. (Error bounds for C 1 -functions). There exists a constant C\ > 0, with C\ < 3 
for Q = S 2 and d < 2 for Q = W% , such that, for f G C 1 ^) : 

||/-Q/||n<Cihw(£)/,V2). 

Proof. Let be the best approximation polynomial of / in P x on the domain O. Consider 
some closed triangular cell T of 7^, n in which we have 

ll/-?*llT = ||/-9*||n. 

Take a point (£, rj) at the midpoint of the external edge of T and let q\ 6 Pi be the linear 
Taylor polynomial of / at that point : 

(8) gi (x, y) = /(£, r?) + P^/(£, r?)(x - + P^/(£, r?)(y - r?). 
Then there hold the following inequalities 

11/ - O/lln < (1 + HQIDII/ - tf-||n = (1 + 11011)11/ - 5*||r < (1 + 11011)11/ - <&||t. 
By Taylor's formula, we have f = q\+r\, with 

(9) n = [P^/M - P^/^K* - + [ D<f>Vf(u,v) - D(°Mf(S,r,)](y - r,), 
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the point (u,v) lying somewhere in the segment joining (£,77) to (x,y). From that, we 
deduce the following upper bound 

\\f- qi \\ T <^(Df,±). 

Actually, for the sake of simplicity, we can assume that T is the triangle with vertices 
(— ?|, 0), (0, |), 0), the point (£,77) being then at the origin. Due to the symmetry of 
the problem w.r.t. the y-axis, we can also assume that (x, y) satisfies x > 0, y > and 
x + y < |. Therefore, as the distance between (u, v) and the origin is bounded above by |, 
we can write : 

|n(x, y)| < W (£>/, ^)(x + y)< ^(Df, ^). 
Finally, as HS2II < 5 and ||W||| < 3, we obtain 

11/ - 5 2 /|| n < 3hu(Df, ^), ||/ - W?/|| n < 2hu{Df, £). 

■ 

Theorem 3. (%) (Error bounds for C 2 -functions). There exists a constant C 2 > 0, with 
C 2 <\forQ = S 2 and C 2 <\forQ = W 2 *, such that, for f G C 2 (Q) : 

\\f-Qf\\n<C 2 h 2 u;(D 2 f, h/2). 



(ii) (Error bounds for C 3 -functions). There exists a constant C 3 > 0, with C 3 < | /or 
Q = S 2 and C 3 <± forQ = , such that, for f G C 3 (n) : 

\\f-Qf\\n<C 3 h 3 \\D 3 f\\. 

Proof. By using a similar technique as in the proof of theorem 2, bounds on Taylor re- 
mainders can be obtained for orders 2 and 3. We have respectively, for / G C 2 (Q) and 
/ G C 3 (tt) 

(10) f = q 2 +r 2 and f = q 3 + r 3 , 
with 

(11) r2 = ^(j P a /(«> «) - y)l 0* - Q1 (y - , 

(12) rs = ^ E QW/K^-o^-^r- 

M=3 

Puting (£, 77) at the origin gives 

1/ - = h < ^v, |) e C) = r ^ 2/ > + ^ ± y^ 2/ > ^ 

|a|=2 V / - 
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/II, 



from which we deduce respectively 

\\f-S 2 f\\<^u(D 2 f,-) and ||/ — W 2 f\\ < —io(D 2 f, ^). 
For / G C 3 (0), we have : 

|/-gs| = |r 3 | < l\\D 3 f\\ PV 1 ^ 2 = l\\D 3 f\\(x + y) 3 < ^\\D 

\a\=3 VO/ 

and finally we obtain 

||/-5 2 /||<^||£> 3 /ll and ||/-^/||<^||D 3 /||. 

(Notice that the inequalities obtained for W 2 are better than those given in Chui-He [2] for 
W 2 ). ■ 

Remark. The constants of the error bounds obtained for S 2 are greater than the corre- 
sponding ones for W%- These results do not mean that S 2 is worse than W 2 , because they 
are a consequence of the fact that S 2 and W 2 belong to the same spline space, i.e. they are 
defined on the same triangulation T mn . Therefore N$ < Nyy, in particular from © and JSJ) 
we have that N$ = 0(mn) and = 0{2mn). 

Now if we assume that S 2 is defined on the triangulation T mn and that W 2 is defined on 
another triangulation obtained by a decomposition of Q into \-^\ ■ subrectangles, 
then the numbers of data values requested by both QIs are almost equal and the constants 
appearing in the respective error bounds are also comparable. 

5 Error analysis for partial derivatives 

In this section, we compute error bounds for the first partial derivatives (Subection 5.2) of 
the quasi-interpolant Qf in S7, and of its second partial derivatives (Subection 5.3) in the 
interior of each triangular cell Ti of T mn . 

5.1 Technical lemmas 

Lemma 1. Let T be a triangular cell of T mn included in the rectangular cell O rs centered 
at M rs , then: 



' 4(/i r )-°i (k s )~ a2 , for |a| = 1 and (x,y) G T 



(13) |£> a %(*,2/)l < < 

(i,j)eK(T) [ 6(h r )- ai (k s )~ a2 , for |a| = 2 and (x,y) G int(T) 

Proof. In the case \a\ = 1, since D a Bij is a linear polynomial in the triangle T = ABC, we 
have: 

\rPBij(x,y)\ KwBxQlPBijiA)^ \D a By(B)\, \D a Bij(C)\}. 
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For \a\ = 2, then D a Bij is a constant inside T. 

In ^3] , the values of the first partial derivatives of B^j at the vertices of T mn and the values 
of the second partial derivatives of Bij inside each triangle T of their support, have been 
computed. Using those values, we can easily deduce the inequalities (fT3*|) . I 

Lemma 2. Let Q be the spline operator given by Let r\ and r2 be the expressions 
defined by and Ml}) for f £ C 1 (r2) and f £ C 2 (ft), respectively. Then for every triangle 
T of T mn , the following majorations hold : 

max <C{hu(Df,h/2), 

(i,J)£K(T) 

where C[ < 30 for Q = S 2 and C[ < 35/2 for Q = W|; 

max |Ai,-(r 2 )| < C 2 h 2 uj(D 2 f, h/2), 

where C' 2 < 61/2 for Q = S 2 and C' 2 < 65/4 for Q = W%. 

Proof. We prove the desired results in the case of the triangle T = T3 (fig. 1) and we use the 
notations of the proof of theorem 1. For the other three types of triangles of T mn we obtain 
the same results using a similar proof, therefore we don't report here the corresponding 
computations. If Q = S 2 , then we know that 

l/%( r i)l = \bijri(Mij) + ainiMi^j) + Cjri(Mj + ij) + %n(Mij_i) + Cjri(Mi ij+1 )\ 

< 3\ ri (M id )\ + |[|n(Afi_i^)| + |n(Afi+ij)| + |n(Afij-_i)| + \ ri (M itj+1 )\ . 

Moreover, taking the origin at the midpoint (Co!??o) of the lower edge of T, we can write 

ri (M id ) = [D^f(Mij) - D^f(0)]( Si - £ ) + [D^f(M id ) - D<M f (0)](tj - r, ), 

where Mi a is some point lying in the segment OMij. 
Finally from the first column of Table 1, we can write 

max |/%(ri)| =max| |/i r , s ±i(ri)|, |/i r± i s (n)|, |// r s (ri)|, |p r± i s _i(ri)| [ = max \n k {r{)\. 

Recall that each central point M k has four neighbours denoted respectively Sk,Ek, Wk- 
Here the biggest constant is obtained for k = 1 corresponding to the central point M rjS+ i. 
In that case, we have 

Mn) < 3|r x (Mi)| + |(|ri(iVi)| + |n(5i)| + |ri(^i)| + |ri(Wi)|), 

where, as ||OMfc|| < ||OMfc||, we can write for example 

„,„,\.,„, 3/i ,„ „ h, 9h ,„ „ h s 



\ r l 



(Afx)| < -D^f(0)\ < -u(Df,3-) < -u{Df,-) 
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|n(Wi)| < h\D^f(W 1 ) - D^f(0)\ + ^{D^fiWx) - D^f(0)\ 

\rx{W 1 )\ < ^u(Df,Vl3^) < 10hu(Df,l). 

From these inequalities and similar ones associated with the three other neighbours of Mi, 
we obtain : 

, /27 1 /25 1\\ , , , h, 

\m(ri)\ < i — + - i- + 20 + - J J fca/p/, -) = 30hu(Df, -). 

In a similar way, we can obtain 

^ 2 (n) and ^ 4 (n) < hu(Df, -), ^(n) < y hu(Df, 

109 /i /i 

/x 5 (ri) and /x 7 (n) < — hu(Df, -), /x 6 (n) < hu{Df,-). 

Finally, we obtain 

max _|Atfe(ri)| < 30huj(Df,-). 
If Q = W?, then from (J3J) and © we get 



|/x*.(n)| < 2| n (M*.)| + T[|n(4*-i.,-i)l + + ki^-r)! + \ri(A*. 



By a procedure similar to that adopted for S 2 , we can obtain that 

max .1^(^)1 <-hu(Df,h/2). 

(i,j)<EK(T) 2 

Now we consider Ay(r2). If Q = S 2 , we know that 

|Wi(r 2 )| < 3|r 2 (M ii )| + -{\r 2 {M^i yj )\ + \r 2 (M i+hj )\ + \r 2 {M iyj _i)\ + \r 2 (M iJ+1 )\), 



where 



r 2 (M M ) = l -{[D^f{M hJ ) - D^f(0)]( Si - £o) 2 + [D^f(M hJ ) - D^fO)]^ - Vo ) 2 

+2[D^f(M id ) - D^f(0)]( Si - £ )(tj - m)h 
with Mij S OMij. Then using a scheme similar to that proposed for r%, we get 

61 

, ma * l/%( r 2)| < —h 2 u{D 2 f,h). 

(l,j)eK(T) z 
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Similarly if Q = W 2 we can deduce: 

M^)) + |r 2 (A* .)|]} < f /^pV, V2). ■ 

Lemma 3. Let Q be a spline operator defined by Q) and let r 3 be defined by Mty) for 
f £ C 3 (ft). Then for every triangle T of T mn , one has the following majoration : 

(14) max r |A^(r 3 )| <C 2 ^\\D^f\\ 
with C' 2 ' < 269/48 for Q = S 2 , C' 2 ' < 65/24 for Q = W£ . 

Proof. As in the proof of Lemma 2, we only discuss here the case of a triangle T = T 3 and 
we don't report the computations for the three other types of triangles. 
If Q = S 2 then, from (|12j) . we can write 

(15) r 3 (M y ) = 1 Yl ( 3 W/(M M -)( Sl -Co) ai (h -r?o) Q2 , 

|«|=3 VQ?/ 

with Mi j S OMij. Therefore if we denote 

(16) p(*,<i) = i E ( 3 )(^-^r(^-%) 

1-1=3 w 
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and if we proceed as in the proof of Lemma 2, from 0, (|15j) . (|16j) we can deduce the desired 
result. Indeed: 

_ max \mj(r 3 )\ < ±||L> 3 /||{39j(,s r ,i s+ i) + ±[ip{s r -i,t s ) + y>(s 
(i,j)€K{T) 

+<p(s r ,t s+2 ) + tp(8 r ,t a )]} < §A 3 ||Z> 3 /IK£[3(i) 3 + i] + 3(f) 3 } < if A 3 ||I> 3 /||- 
Using the same method, we can prove (|14j) for Q = W 2 *.B 

5.2 Error estimates on first partial derivatives 

Theorem 4. (C 1 functions). For \a\ = 1, there exists a constant C\ > 0, with C\ < 120 
for Q = S 2 and C x < 70 for Q = W| , sucfc toaf, /or / G C 1 ^) 

(17) ||^ a /-^ a Q/lb< [l+Ui(^)]a;(I?/,V2) 

Proof. For |a| = 1 we consider a closed triangular cell T of T mn where 

\\D a f - D a Qf\\ n = \\D a f - D a Qf\\ T . 
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For any point P = (x,y) £ T, since Q reproduces P2, we can write: 

(18) \D a f(P) - D a Qf(P)\ < \D a f(P) - D a qi (P)\ + \D a Q(f - qi )(P)\, 

with q± denned in jSJ). 

We remark that, from (JHJ), there results: 

(19) \D a f(P) - D a gi (P)\ = \D a f(P) - D a f{0)\ < u(Df, h/2) 
and 

\D a Q(f - qi )(P)\ < max |Ao-(ri)| £ ^B^P)] 
with r\ defined in @. 

Moreover we recall that, from Lemma 1, for |a| = 1 

(20) Yl l^ a ^i(^)l <4r 2 . 

(<d')e#Cr) 

Finally from Lemma 2 and (|2U|) we obtain 

(21) \D a Q(f- gi )\ <^u(Df,h/2) 

Therefore, from (P |) . (fTD |) and ipTJl . the result JTZJ) follows with C x = AC[. I. 

Theorem 5. (i) (C 2 functions). There exists a constant C 2 > 0, with C 2 < 122 for Q = S 2 
and C 2 < 65 for Q = WJ, such that, for f G C 2 (Q) and \a\ = 1 : 



(22) ||I> a /-£) Q Q/l|n< 



1 + C 2 



hcv(D 2 f, h/2) 



(ii) (C 3 functions). T/ien there exists a constant C3 > 0, with C3 < ^ /or Q = S 2 and 
C 3 <f forQ = W 2 *, suc/i t/iat, /or / G C 3 (ft) and |a| = 1 : 

(23) ||Zr/-D a Q/||<C 3 (j)/i 2 ||^/ll 

Proof. The proof is similar to that of Theorem 4. 
For |a| = 1 and P £ T, we can write: 

(24) |£> Q /(P) - ^ Q Q/(P)| < |D°/CP) - ^ 2 (P)| + \D a Q{f - g 2 )(P)| 
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where q 2 G P2 has been defined in (fTU)) . 
We remark that 

(25) \D a f(P) - D a q 2 (P)\ < hcv(D 2 f, h/2). 



Moreover 

(26) \D a Q(f - q 2 )(P)\ < max . |Ay(r 2 )| V |D Q ^(P)|, 

(ij)£i " (T) (i,i)^(T) 



with r 2 defined in IJ11|) . 

Now from l)26|) . Lemma 1 and Lemma 2 we can write that 



(27) 



\D a Q(f-q 2 )(P)\<4C' 2 (jj huj(D 2 f, h/2). 



Finally from (21, J2SJ) and (g7|), the result Q follows, with C 2 = 4C 2 . 
Moreover if / G C 3 (J1), from (|12|) . there results: 



(28) 



\D a (f-q 3 )(P)\<^-\\D 3 f\\, 



and from Lemma 1 and Lemma 3 

(29) \D«Q(f - q 3 )(P)\ < AC 2 ' (jj h 2 \\D 3 f\\. 

Therefore from dHJ), (gHD and © the result (22J) follows, with C 3 = 4C 2 '. 
5.3 Error estimates on second partial derivatives 



Theorem 6. (i) (C functions). There exists a constant D 2 > 0, with D 2 < 183 for 
195 t,,.. r\ _ ri-i ..„..;, //.,,/ r.,,, r / ''-( 



Q = 5 2 , and D 2 < ^ for Q = W$, such that for f G C 2 (Q) and |a| = 2 ; 

2. 



(30) 



|^ a /-^Q/||int(T) < 



1 + Do I % 



to(D 2 f,h/2). 



(ii) (C 3 functions). There exists a constant D% > 0, with D3 < if Q = S 2 and D% < 



65 



i/Q = W$, such that, for f G C 3 (0) : 

(31) \\D a f-D a Qf\U {T) < 



I + D3 



Proof. For |a| = 2 and for any P G int(T), we have 

(32) \D a f - D a Qf\ < \D a f - D a q 2 \ + \D a Q(f - q 2 ) 

with q 2 defined in (I1UI) , 
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Prom ifTUjl . we deduce 

(33) \D a f-D a q 2 \<io(D 2 f,h/2). 
Moreover 

(34) \D a Q(f-q 2 )\ < max |Ay(r 2 )| ■ V |£> a ^|, 

{JJ) ^ (T) d,j)eK(T) 

with r2 defined in ((lip. 

From (|34|). Lemma 1 and Lemma 2 we obtain that for P eint(T): 

(35) \D a Q(f-q 2 )(P)\ <6C 2 (J\ u(D 2 f,h/2) 

Therefore from (|32J, © and © the result dSUJ follows, with D 2 = 6C' 2 . 

Finally if / € C 3 (fi), for P G int(T) and |a| = 2, from (|32|) . Lemma 1 and Lemma 3, we 

obtain the result (f3"TJ). with D 3 = 6C 2 '. I 

We note that the same remarks given at the end of Section 3 are also valid for the error 
bound constants of the above theorems. 

5.4 Convergence for quasi- uniform partitions 

Assume that the sequence of partitions {X m x Y n } of Q is 7-quasi uniform i.e. there exists 
a constant 7 > 1 such that 

< h mn / 5 mn ^ 7) 

where /i mn and 5 mn are respectively the maximum and the minimum steplengths of the 
partition {X m x Y n }. Then the following theorem shows that for both Q = S 2 and Q = W 2 

D a Qf^D a f as h mn ^0 

in fl for I a I = 1, and in the interior of each triangular cell T of T mn , for \a\ = 2. 

Theorem 7. Let {X m x Y^} a "f-quasi uniform sequence of partitions. 

(i) If f € C s (^), a = 1,2, t/ien /or |a| = 1 

HIT/ - D a Qf\ \q = 0(h s m MD s f,h mn /2)) 

(ii) /// € C 2 (ft) then for\a\ = 2 

||D<7 _ D a Qf\\ int(T) = 0(u;(D 2 f, h mn /2)) 
Proof. The result immediately follows from the 7-quasi uniformity of {X m x Y n } and 
from Theorems 4, 5, 6. I 
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